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We propose two types of topologically stable knot solitons in condensed matters, one in two- 
component Bose-Einstein condensates and one in two-gap superconductors. We identify the knot in 
Bose- Einstein condensates as a twisted vorticity flux ring and the knot in two-gap superconductors 
as a twisted magnetic flux ring. In both cases we show that there is a remarkable interplay between 
topology and dynamics which transforms the topologcal stability to the dynamical stability, and vise 
versa. We discuss how these knots can be constructed in the spin-1/2 condensate of 87 Rb atoms 
and in two-gap superconductor of MgBi. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Vs, 74.72,-h 
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I. INTRODUCTION 

Topological objects (monopoles, vortices, skyrmions, 
etc.) have played an increasingly important role in 
physics Q, 0- In particular finite energy topological 
objects have been widely studied in theoretical physics 
0. LI A recent addition to this family of finite energy 
solitons has been the knots @, |(| . The interest on these 
topological objects, however, has been confined mainly 
to theoretical physics, because most of them (exept the 
vortices) exist in "hypothetical" worlds which are very 
difficult to create in laboratories. The only topological 
objects which one can realistically expect to exist in the 
"standard" models are the electroweak monopoles and 
dyons in Weinberg-Salam model which have a non-trivial 
and Z M dressing 0,0- Unfortunately these objects 
could carry an infinite energy, which makes it impossible 
to create them experimentally. 

Fortunately the recent experimental advances in con- 
densed matter physics, in particular the construction 
of multi-component Bose-Einstein condensates (BEC) 
made of 87 Rb |^.ITol| and two-gap superconductors made 
of MgBi [Til ll^ , have widely opened new opportunities 
for us to create new topological objects experimentally 
which so far have been only of theoretical interest. The 
purpose of this paper is to argue that these new multi- 
component condensates could allow us to have real knots, 
topologically stable finite energy 3-dimensional solitons. 
This is because, due to the multi-component structure, 
the new condensates have a non-Abelian symmetry which 
provides the needed topology for the stable knots. 

To understand how the realistic knots can appear in 
these condensed matters, it is necessary to understand 
the prototype Faddeev-Niemi knot in Skyrme theory. 
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The Skyrme theory is well known to have a magnetic vor- 
tex known as the baby skyrmion [l3j . and the Faddeev- 
Niemi knot can be identified as a twisted vortex ring 
made of a helical baby skyrmion [y, Il4| . In the following 
we show how one can construct similar knots from heli- 
cal vortices in two-component Bose-Einstein condensates 
and two-gap superconductors. 

The paper is organized as follows. In Section II we 
discuss how the helical vortex can give rise to the proto- 
type knot in Skyrme theory for later purpose. In Section 
III we discuss gauge theory of two-component BEC which 
has the vorticity interaction, and show that the theory is 
closely related to Skyrme theory. With this we show that 
the theory allows a topological knot similar to the one in 
Skyrme theory. We also identify that this knot is a vortic- 
ity knot very similar to the one in Gross-Pitaevskii the- 
ory of two-component BEC. In Section IV we discuss the 
Landau-Ginzburg theory of two-gap superconductor and 
show that the theory is closely related to the gauge the- 
ory of two-component BEC. With this we argue that the 
theory can also admit a knot, a twisted magnetic vortex 
ring, similar to the vorticity knot in two-component BEC. 
In doing so we also establish the non-Abelian flux quan- 
tization in two-gap superconductor. We demonstrate the 
existence of magnetic vortex whose flux is quantized in 
the unit 4-ir/g, not 2ir/g, in two-gap superconductor. Fi- 
nally in Section V we discuss physical implications of our 
results. 



II. KNOT IN SKYRME THEORY 



The Skyrme theory has a rich topological structure. 
It has skyrmion and baby skyrmion 0, El- But re- 
cently it has been shown that the theory allows a knotted 
soliton identical to the Faddeev-Niemi knot in Skyrme- 
Faddeev non-linear sigma model, which can be identified 
as a twisted vortex ring made of helical baby skyrmion 
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[^.Il4|. To see this, we review the knot in Skyrme theory 
first. 

Let ui and n be the scalar field and the non-linear 
sigma field in Skyrme theory. With 



U = exp(u; — • n) = cos — — iy<j ■ n) sm — , 
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the Skyrme Lagrangian is expressed as 
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£ = ^tr /,, 2 , -^tr( / ./ 
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-(d^uj) 2 + {l-cosu;)(d l ,h) 2 

1 - COSUJ 2 

(d^ujd^n ~ OvUO^n) 



+ (1 — cosw) 2 (<9 M n x dyfif 
The equation of motion is given by 
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+ — (1 - cos uj)d 'J (d^ojd^n - d v ujd^h) ■ d u h] 
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(1 — coscl>) smuj(d^n x d v fi) 2 = 0, 
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— cosw)ti x d^n 

+ ^(1 — cos uj) [(d u uj) 2 n x dfj,n — (d^ujd v uj)h x d u n\ 

+ ^(1 — cosw) 2 (n • d^n x d v h)d u fi\ = 0. (3) 
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With the spherically symmetric ansatz 

u> = uj(r), h = f, (4) 

(3) is reduced to 

rsin 2 (w/2) d 2 uj 



cPoj 2 dio 2 sin ui a 
dr 2 r dr r 2 u 2 



sinw duj 2 2sinwsin 2 (o;/2) 



2r 2 ^dr^ 



dr 2 
0. 
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FIG. 1: The baby skyrmion (dashed line) with m — 0,n — 1 
and the helical baby skyrmion (solid line) with m = n = 1 in 
Skyrme theory. Here q is in the unit y/a/ p and k = 0.8 y^a. 



Notice that since forms a closed two-form, it always 
admits a t/(l) potential C M . 

The equation Q allows non-Abelian monopole, baby 
skyrmion, helical baby sky rmion, and Faddeev-Niemi 
knot as its solutions |g, [Tj| . But for our purpose it is 
important to understand the helical baby skyrmion, be- 
cause this plays a crucial role for us to construct the knot. 
So we review the helical baby skyrmion. 

To construct the desired helical baby skyrmion let 
(g, ip, z) the cylindrical coodinates, and choose the ansatz 



' sin f(g) cos (mkz + rup) \ 
sin f(g) sin (mkz + nip) \ 
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cos f(g) I 
With this the equation @ is reduced to 

/ / t , t n 2 , sin 2 f \ ■■ / 1 p 
1 + (m 2 k 2 + —)^-L)f + - + 2 1 - 
\ g z g p ' \g p 

. 9 , 9 n 2 , sin / cos / ; 1 , 9 , 9 n 2 , sin 2 f \ ■ 
+(m 2 k 2 + —) i^±f +-{m 2 k 2 - — )- r L / 
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i sin / cos / = 0. 



Imposing the boundary condition w(0) — 2tt and w(oo) = So with the boundary condition 
0, one has the well-known skyrmion 



A remarkable point of Skyrme theory is that ui — tt, 
independent of h, becomes a classical solution So 
restricting u to tt, one can reduce the Lagrangian to 
the Skyrme-Faddeev Lagrangian 



/(0)=7T, /M=0, 



(9) 



(10) 



C -> -^-(a M n) 2 - T (9 M n x a„fi) 2 
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in which case the equation of motion @ is reduced to 



n x <9 2 n + —{d ll N lllJ )d u h = 0, 
A* 

N MV = h ■ (d^n x d v n) = ^C,, - d^C^. 



(7) 



we obtain the non-Abelian vortex solutions shown in 
Fig^ Notice that, when m = 0, the solution describes 
the well-known baby skyrmion [l3j. But when m is not 
zero, it describes a helical baby skyrmion, a twisted mag- 
netic vortex which is periodic in z-coodinate |14| . In this 
case, the vortex has a quantized magnetic flux not only 
along the z-axis but also around the z-axis. 

The helical baby skyrmion will become unstable un- 
less the periodicity condition is enforced by hand. But 
it plays a very important role because one can make it a 
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vortex ring by smoothly connecting two periodic ends. In 
this case the vortex ring acquires the topology of a knot, 
and thus becomes a knot itself 00] . In fact it becomes a 
knot made of two magnetic fluxes linked together, whose 
knot topology is described by the Chern-Simon index of 
the potential C M , 



Q = 



32tt 2 



1k CiN lk d 3 x = ran, 



(11) 



which describes the Hopf mapping ir^S 2 ) defined by h. 

Clearly the knot has a topological stability, because 
two flux rings linked together can not be disjointed by a 
smooth deformation of the field configuration. Moreover 
the topological stability is backed up by a dynamical sta- 
bility. This is because the knot can be viewed as two 
magnetic fluxes linked together, and the magnetic flux 
trapped in the knot disk can not be squeezed out. This 
provides a stablizing repulsive force which prevent the 
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collapse of the knot 

From our discussion it becomes clear that the ex- 
istence of the helical baby skyrmion is crucial for the 
existence of a topologically stable knot. In the fol- 
lowing we show that an identical mechanism works in 
two-component BEC and two-gap superconductor which 
guarantees the existence of a stable knot. 



III. KNOT IN TWO-COMPONENT 
BOSE-EINSTEIN CONDENSATES 



The recent advent of multi-component BEC (in par- 
ticular the spin- 1/2 condensate of 87 Rb atoms) has 
widely opened a new opportunity for us to study novel 
topological objects which can not be realized in ordinary 
(one-component) BEC 0,0] . This is because the multi- 
component BEC naturally allows a non-Abelian struc- 
ture which accommodates a non-trivial topological ob- 
jects, in particular a topolgical knot which is very similar 
to the knot in Skyrme theory [TH Hi} . 

There are two competing theories of two-component 
BEC, the popular Gross-Pitaevskii theory [T^j and the 
gau ge theory of two-component BEC proposed recently 
Both theories predict topological knots. Many au- 
thors have already claimed the existence of a knot in 
Gross-Pitaevskii theory [T^j . Moreover it has been shown 
that this knot can be identified as a vorticity knot which 
is made of two vorticity fluxes linked together, whose 
topology %3(S 2 ) is fixed by the Chern-Simon index of 
the velocity potential of the condensate |l6( • 

So in this section we discuss the gauge theory of two- 
component BEC, and show that this theory also has a 
vorticity knot very similar to the knot in Gross-Pitaevskii 
theory. Consider a "charged" two-component BEC de- 
scribed by a complex doublet <f> interacting "electromag- 
netically" , which can be described by the following non- 



relativistic gauged Gross-Pitaevskii type Lagrangian (we 
will discuss the relativistic generalization later) 



tf(D t <t>) - {Dt4>)U 
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1 -f 2 
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where D^<f> = (d^ — igA^tfi, and fj? and A are the chemi- 
cal potential and the quartic coupling constant. Normal- 
izing (j> to (\/2m/h)(f> and putting 



= 1) 



(13) 



we have the following Hamiltonian in the static limit, 



H 



1 



(d lP f 



V 2 2 



8 P 



If 2 

4 y 



(14) 



Notice that here we have also normalized /i 2 and A to 
(h 2 /2m)/Lt 2 and (H 2 /2m) 2 A. From now on we will use the 
normalized Hamiltonian (11411 . 

At this point it is important to realize that the "elec- 
tromagnetic" interaction here should be self-induced, 
since we are dealing with neutral condensates. So we 
identify the "electromagnetic" potential by the velocity 
field of the doublet £, 



4 = 

9 



(15) 



A justification for this is that the U(l) gauge invariance 
almost dictates us to identify the velocity field as the 
gauge potential. Indeed in the absence of the Maxwell 
term, (| 1 o P becomes an equation of motion. With this 
identification the field strength becomes the vorticity 
field 



- fj v 

i 
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(16) 



One might wonder why we have included the vortic- 
ity interaction in the Lagrangian (|12Jl . when this is ab- 
sent in the Gross-Pitaevskii Lagrangian. The reason is 
because it costs energy to create the vorticity in two- 
component BEC. In ordinary BEC one does not have to 
worry about the vorticity because the vorticity is iden- 
tically zero (since the velocity is given by the gradient 
of the phase of the one-component complex condensate) . 
But a non-Abelian (multi-component) BEC has a non- 
vanishing vorticity, in which case it is natural to include 
the vorticity interaction in the Lagrangian [TEl ITtij . 

With HI 5(1 the Hamiltonian ((14(1 becomes 



H 



(M 2 + |-(l^| 2 -|^| 2 )-|(p 2 -^) 2 
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Po = 



2jS 

T 



(17) 



Notice that now the coupling constant g represents the 
coupling strength of the vorticity. Minimizing the Hamil- 
tonian we have the following equation of motion 



1 



g 2 p 2 

= 0. (18) 
To understand the meaning of this notice that with 



we have 



h = 



\dM 2 -\^d,e = 7 (d,nf 



i{d^d u £ - d v ^d^) = -h ■ (d^ii x d^fi) 



(19) 



(20) 



where N^ v is mathematically identical to what we have 
in Skyrme theory in Q. This tells that the potential 
for the two-form in @ is given by (up to a gauge 
transformation) 

C M = 2gA il = -2i^d^. (21) 

More importantly with H20|l the Hamiltonian l|17l) can be 

expressed as 

,>2 



n= 1 -(d l pf + ^(d z n) 2 + ^p 2 



Pi? 



1 



■{difi x djfi) 



(22) 



So the theory becomes very similar to Skyrme-Faddeev 
theory, which strongly indicates that the gauge theory 
of two-component BEC can allow a knot. As impor- 
tantly this strongly implies that the Skyrme-Faddeev the- 
ory could play an important role in condensed matter 
physics. 

To simplify the equation (|18fl notice that from i|15|) 
and l|19fl we have 



= (igAfi + —(7 ■ d^h)^. 



(23) 



With the identity the second equation of H18|) is reduced 
to 



B = d 2 n + 



(A + B ■ <r)£ = 
A-(5 4 n) 2 , 

i 

-c) l n - — 

P 



2g 2 p 



-^diNijdjh, 



(24) 



which can be written as 

A + B-h = Q, 
h x B — in x (h x B) = 0, 



(25) 



n x (d 2 h + 2—d.n) + — — ^jV^^n = 0. (26) 
P ^9 2 P 2 

So we can put l|18f) into the form 

d 2 p-\{d % n) 2 p=\{p 2 -pl) Pl 



x d 2 



2 n x dj 



1 



g 2 p 2 



diN i3 d 3 h = Q. (27) 



This is the equation of two-component BEC that we are 
looking for. The similarity between this and the equation 
Q) of Skyrme-Faddeev theory is unmistakable. 

Notice that the target space of £ and n is S 3 and S 2 , 
but here we have transformed the equation for £ in l|18|) 
completely into the equation for h in (|27|l . This is made 
possible because, with the Abelian gauge invariance of 
(6), the physical target space of £ becomes the gauge 
orbit space S 2 — S^/S 1 which forms a CP 1 space. This 
means that we can view the theory as a self interacting 
CP 1 theory (coupled to a scalar field p), and replace £ 
completely in terms of n. 

To show that the theory has a knot solution we con- 
struct a helical vortex solution first. To do this we choose 
the ansatz 



p = p(q), £ 



Ho) , ■ v 
cos — - — exp(— inip) 



sin 



2 

fie) 



exp(imkz) 



h = fto£ = 



' sin f(g) cos (mkz + nip) \ 
sin f(g) sin (mkz + nip) . 
cos /(p) / 



(28) 



With this (H3 is reduced to 

1 . 1 

Q 



P+-p-i[f 2 + (m 2 k 2 + ^ T )sm 2 f)p 



\{p 2 - pDp, 



, „ , „ n 2 . sin 2 f 
l + (m 2 k 2 + —)- T f 



, ,.,,/+ (± + 2* 
Q 9 P ' v £> P 



n l sin / cos / ; 1. 9l9 n 2 .sin/\„ 



9 P 



g p 



-(m 2 k 2 



:) sin /cos / = 0. 



So with the boundary condition 

p(0) = 0, p(oo) = po, 
/(0) = tt, /(oo) = 0, 



(29) 



(30) 
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This has two vorticity fluxes, 4>z along the z-axis 
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FIG. 2: The non-Abelian vortex (dashed line) with m = 
0, n = 1 and the helical vortex (solid line) with m = n = 1 in 
the gauge theory of two-component BEC. Here we have put 
\/g 2 = 1, k = 0.64 V\p , and q is in the unit of 1/V^Xpo. 



we obtain the non-Abelian vortex solutions shown in 

There are three points that have to be emphasized 
here. First, when m = 0, the solution describes the 
untwisted non-Abelian vortex |15|. But when m is not 
zero, it describes a helical vortex which is periodic in z- 
coodinate. In this case, the vortex has a non-vanishing 
velocity current (not only around the vortex but also) 
along the z-axis. Secondly, the doublet 4> starts from the 
second component at the core, but the first component 
takes over completely at the infinity. This is due to the 
boundary condition /(0) = 7r and /(oo) = 0, which as- 
sures that our solution describes a genuine non-Abelian 
vortex. Thirdly, when / = (or / = it) the doublet 
effectively becomes a singlet, and (|29f) describes the well- 
known vortex in single-component BEC. Only when / 
has a non-trivial profile, we have a non-Abelian vortex. 

In Skyrme theory the helical vortex is interpreted as 
a twisted magnetic vortex whose flux is quantized 0, Ih4| 
Now we show that the above vortex is a twisted vorticity 
flux which is also quantized. To see this notice that the 
non-Abelian structure of the vortex is represented by the 
doublet £. Moreover, the velocity field of the doublet is 
given by 



-4„ 



25 



(cos / + Vjd^ip 



mk 



(cos / - l)<9 M z, 



(31) 



which generates the vorticity 



H liV = --{d^d v t-d v £}d ll 

= ]_ 

2g 



sin f (n(d,j,Qd v ip - d u gd^ip) 



+mk(d fl Qd„z - d v Qd ll z) S j . 



(32) 



H^gdgdip 



2nn 
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(33) 



and <p0 around the the z-axis (in one period section from 
z = to z = 2n/k) 



IT A A 2lTm 

Hzgdgaz = . 

D " 9 



(34) 



This shows that the helical vortex is made of two quan- 
tized vorticity fluxes, the (/>z-flux which is concentrated 
at the core and the </>^-flux which surrounds it. This 
confirms that the helical vortex is a twisted vorticity 
flux which is very similar to the helical vortex in Gross- 
Pitaevskii theory [T^ . 

Now, with the helical vortex, one can easily make 
a twisted vortex ring smoothly connecting the periodic 
ends together. And just as in Skyrme theory the twisted 
vortex ring becomes a topological knot. But here it is 
TT^iS 3 ) of the doublet £ which provides the non-trivial 
quantum number, 



1 

472 



iktfdi&djgdkSd 3 



(35) 



Of course this is identical to the expression Qllfl. due to 
the Hopf fibering of S 3 to S 2 x S* 1 16] . This tells that we 
can express the knot quantum number either by tt^(S 3 ) 
or by n 3 (S 2 ). 

Obviously this knot has a topological stability, be- 
cause the knot topology can not be changed by a 
smooth deformation of the field configuration. Moreover 
it has a dynamical stability. To understand this notice 
that the knot has a twisted velocity field so that it has 
a non-vanishing velocity around the z-axis. This means 
that it carries a non- vanishing angular momentum along 
the z-axis. And this angular momentum provides the 
dynamical stability, because it creates a centrifugal force 
that prevents the collapse of the knot. Notice that this 
dynamical stability originates from the knot topology, 
because the angular momentum comes from the twisted 
velocity field. In this sense the topological stability and 
the dynamical stability have one and the same origin. 
It is this remarkable interplay between topology and dy- 
namics which assures the stability of the knot. The non- 
trivial twisted topology of the knot expresses itself in the 
form of the angular momentum, which in turn provides 
the dynamical stability of the knot. This presence of the 
angular momnetum is what differentiates our knot from 
the untwisted Abrikosov-type vortex ring which has no 
dynamical stability. 

There have been assertions that two-component BEC 
admits a knot ^(J ■ But notice that this knot is based on 
Gross-Pitaevskii theory of two-component BEC, which 
has no vorticity interaction. In contrast our knot is based 
on the gauge theory in which the vorticity interaction 
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plays a crucial role. Nevertheless physically two knots 
are very similar [3. Both can be identified as a vorticity 
knot. This implies that both theories should be taken 
seriously as a theory of two-component BEC. 



To obtain the desired knot we first construct a super- 
conducting helical magnetic vortex. Let 



p = p(e), £ 



f(o) i ■ v 
cos exp(— imp) 

■ ^) C h ^ 

, sin — - — exp(iTOfcz) , 



IV. KNOT IN TWO-GAP 
SUPERCONDUCTORS 



In the above gauge theory of spin-1/2 condensates 
the gauge interaction was a self-induced interaction. But 
when the doublet is charged, the gauge interaction can be 
treated as independent. In this case the theory can de- 
scribe a two-gap superconductor. But even in this case 
the knot topology and thus the knot itself should sur- 
vive. This implies that two-gap superconductor should 
also have a topological knot. 

The knot in two-gap superconductor could be ei- 
ther relativistic or non-relativistic, and appear in both 
Abelian and non-Abelian setting [T3. In this paper we 
will discuss the relativistic knot in the Abelian setting 
(a non-relativistic Gross-Pitaevskii type theory gives an 
identical result). Consider a charged doublet scalar field 
(j> coupled to the real electromagnetic field, 



C = -\D„ 



) 2 — -F 2 



(36) 



The Lagrangian has the equation of motion 

,2 



D 2 d> = A 



EL 

A 



= j u = ig (D v <j>y<f> - <P(D v <j>) . (37) 



Now, with 



£ f £ = l, & = (38) 



we can reduce (|3*Tj) to ^3 



d 2 p (\{d»hf + g 2 (A, + A,f)p = ±{p 2 - p 2 )p, 



n x d 2 h + x d^h H — ^d^F^d^h = 0, 



2 

p ' gp 2 

d fl F fil/ =j f2 =g 2 p 2 {A f , + A f ,), 

A fl = — Cd^, pa = — . 
9 ^ 



(39) 



This is the equation for two-gap superconductor. Notice 
that with A^ = —A^ the first two equations reduce to 
l|27jl. This tells that the gauge theory of two-component 
BEC and the above theory of two-gap superconductor 
are closely related. 



A^ = -(nAi(g)d ll (p + mkA 2 (g)d fl z). 

sin /(e) cos {tup + mkz) " 
sin /(e) sin (nip + mkz) 
cos /(p) 

(cos /(g) + l)d M y> 
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mk 

V 



cos /(e) 



1)^2. 



(40) 



With this we have 



= 9P 



i(Ai- 



cos / + 1 



du<P 



+mk(A 2 
and l|39|l becomes 



cos / — 1 



(41) 
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Q 



cos / + 1 



m 2 k 2 ) sin 2 / 

cos / — 1 



^ 2 k 2 [A 2 - 



t;(p 2 - pDp, 



1 



/ + (±+2^)/-2ft(A 1 

e p v e 



^ 2 k 2 (A 2 



A 2 



A 1 -gV(A 1 
A 2 -g 2 p 2 [A 2 



sin / = 0, 

cos / + 1 
2 

cos/ 



1 



0, 
0. 



(42) 



Now, we impose the following boundary condition for the 
non-Abelian vortices |l7j . 

p(0)=0, p(oo)= P0 , /(0)=7T, /(oo)=0, 

Ai(0) = -1, Ai(oo) = l. (43) 

This need some explanation, because the boundary value 
AxiO) is chosen to be —1, not 0. This is to assure the 
smoothness of p(g) and /(e) at the origin. Only with 
this boundary value they become analytic at the origin. 
At this point one might object the boundary condition, 
because it creates an apparent singularity in the gauge 
potential at the origin. But notice that this singularity 
is an unphysical (coordinate) singularity which can eas- 
ily be removed by a gauge transformation. In fact the 
singularity disappears with the gauge transformation 



■>exp(inip), A^ — > A^ 



1 a 



(44) 
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FIG. 3: The non-Abelian vortex (dashed line) with m = 
0, n = 1 and the helical vortex (solid line) with m = n = 1 
in two-gap superconductor. Here we have put g = 1, A = 2, 
k — po/W, and g is in the unit of 1/po- Notice that A2 has a 
logarithmic singularity at the origin. 



which changes the boundary condition A\(Q) = — 1 and 
Ai(oo) = 1 to Ai(0) = and Ai(oo) = 2. Mathemat- 
ically this boundary condition has a deep origin, which 
has to do with the fact that the Abelian U(l) runs from 
to 2tt, but the S 1 fiber of SU(2) runs from to 4tt 
[l7j . As for ^2(0), we choose ^2(00) = to make the 
supercurrent vanishing at infinity and require the vortex 
superconducting. As we will see, this requires a logarith- 
mic divergence for ^2(0). The boundary condition will 
have an important consequence in the following. 

With the boundary condition we can integrate l(42|) 
and obtain the non-Abelian vortex solution of the two- 
gap superconductor, which is shown in Fig|SJ The solu- 
tion is very similar to the one we have in two-component 
BEC. When m — 0, the solution (with A 2 = 0) describes 
an untwisted non-Abelian vortex |17j . But when m is not 
zero, it describes a helical magnetic vortex which is pe- 
riodic in z-coordinate. Moreover, the vortex starts from 
the second component at the core, but the first compo- 
nent takes over completely at the infinity. This is due to 
the boundary condition /(0) = n and /(oo) = 0, which 
assures that our solution describes a genuine non-Abelian 
vortex. This is true even when m = 0. Only when / = 
(or / = tt) the doublet effectively becomes a singlet, and 
(|42|l describes the Abelian Abrikosov vortex of one-gap 
superconductor. 

There are important differences between the non- 
Abelian vortex and the Abrikosov vortex. First the non- 
Abelian vortex has a non-Abelian magnetic flux quanti- 
zation . Indeed the quantized magnetic flux </> 2 of the 
non-Abelian vortex along the z-axis is given by 



H z = 



n A x 



I 



H z d 2 x 



9 Q 
2im 



[Ai(oo)-Ai( 



Aim 
9 



(45) 



Notice that the unit of the non-Abelian flux is Air/g, not 
2ir/g. This is a direct consequence of the boundary con- 
dition H43J1 ■ This non-Abelian quantization of magnetic 
flux comes from the non-Abelian topology ir 2 (S 2 ) of the 
doublet £, or equivalently the triplet h, whose topological 
quantum number is given by 



tripl' 

m 



q = h\ Hz<fx = ~~Lj ^Atdj^dPx 



eijfi ■ (difi x djfi)d 2 x = n. 



(46) 



This distinguishes our non-Abelian vortex from the 
Abelian vortex whose topology is fixed by 7r 1 (5' 1 ). 

Another important feature of the non-Abelian vortex 
is that it carries a non-vanishing supercurrent along the 
z-axis, 



i z = mkg J p 2 (A 



cos f + I 



) gdgdip 



2irmk 
9 



9 

{qM) 



(A 2 + -A 2 )gdg 



g—oo 



9 



(47) 



This is due to the logarithmic divergence of A 2 at the 
origin. 

Notice that the superconducting helical vortex has 
only a heuristic value, because one needs an infinite en- 
ergy to create it (since the magnetic flux around the vor- 
tex becomes divergent because of the singularity of A 2 
at the origin). With the helical vortex, however, one can 
make a vortex ring by smoothly bending and connect- 
ing two periodic ends. In the vortex ring the infinite 
magnetic flux of A 2 can be made finite making the finite 
supercurrent 1)47(1 of the vortex ring produce a finite flux, 
and we can fix the flux to have the value AfKmj g by ad- 
justing the current with k. With this the vortex ring now 
becomes a topologically stable knot. 

To see this notice that the doublet £, after forming a 
knot, acquires a non-trivial topology 113 (S 2 ) which pro- 
vides the knot quantum number, 

Q = -4^2 j '..^C'>^''A O^utKr 
"' 1 e ijk C i (d j C k -d k C j )d 3 x = mn. (48) 



32tt 2 



This is nothing but the Chcrn-Simon index of the poten- 
tial C^,, which is mathematically identical to the quan- 
tum number of the knots we discussed before. This tells 
that our knot is also made of two quantized magnetic 
flux rings linked together whose knot quantum number is 
fixed by the linking number mn. Obviously two flux rings 
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FIG. 4: The regular helical magnetic vortex with m = n = 1 
in two-gap superconductor. Here we have put g = 1, A = 2, 
k = 0.12/9 , ^2(0) = 0.5, and g is in the unit of 1/po- Notice 
that the solution is completely regular. 



linked together can not be separated by any continuous 
deformation of the field configuration. This provides the 
topological stability of the knot. 

Again this topological stability is backed up by a dy- 
namical stability. To see this notice that the supercurrent 
of the knot has two components, the one around the knot 
tube which confines the magnetic flux along the knot, but 
more importantly the other along the knot which creates 
a magnetic flux passing through the knot disk. This com- 
ponent of supercurrent along the knot now generates a 
net angular momentum which provides the centrifugal re- 
pulsive force preventing the knot to collapse. This makes 
the knot dynamically stable. 

To compare our knot with the Abrikosov vortex ring 
(made of the Abrikosov vortex in conventional supercon- 
ductor), notice that the Abrikosov knot is empty (i.e., 
does not carry a net supercurrent). As importantly it 
is unstable, and collapses immediately. In contrast our 
knot has a helical supercurrent, and is stable. Further- 
more these two features are deeply related. The heli- 
cal supercurrent plays a crucial role to stablize the vor- 
tex ring by providing the net angular momentum, which 
prevents the collapse of the vortex ring. And this heli- 
cal supercurrent originates from the knot topology. This 
remarkable interplay between topology and dynamics is 
what provides the stability of the knot. The nontrivial 
topology expresses itself in the form of the helical super- 
current, which in turn provides the dynamical stability of 
the knot. We emphasize that this supercurrent is what 



distinguishes our knot from the Abrikosov vortex ring, 
which has neither topological nor dynamical stability. 
V. DISCUSSION 



In this paper we have presented a compelling ar- 
gument for the existence of topological knots in two- 
component BEC and two-gap superconductor. Similar 
knots have popped out almost e very where, in particular 
in high energ y p hysics in QCD (l8j and Weinberg and 
Salam model [lJj- But we emphasize that at the center 
of these topological objects lies the baby skyrmion and 
the Faddeev-Niemi knot 0,^3- I* 1 fact, our helical vor- 
tices and knots in this paper are a straightforward gen- 
eralization of the baby skyrmion and the Faddeev-Niemi 
knot. 

It has been assumed that the topological objects in 
Skyrme theory can only be realized at high energy, at 
the hadronic scale. But our analysis shows that simi- 
lar objects could exist in a completely different environ- 
ment, at a much lower scale, in low energy condensed 
matters. If so, the challenge now is to verify the exis- 
tence of the topological knot experimentally in condensed 
matters. Constructing the knot might not be so easy at 
present moment. Nevertheless, with some experimental 
ingenuity, one should be able to construct the knots in 
condensed matters. 

Note Added: One might doubt the existence of a su- 
perconducting knot because the superconducting helical 
vortex we discussed in Section IV was singular (and thus 
unphysical). In this note we report a regular supercon- 
ducting helical vortex which has a finite magnetic flux 
around the axis and thus a finite energy. The regular so- 
lution is obtained linking ^(0) with k. For example for 
k = 0.12 we obtain the regular solution shown in FigQ] 
with ^2(0) = 0.5. This type of regular helical vortex 
has vanishing supercurrent i z , but could still be called 
superconducting because it has a non-trivial supercur- 
rent density j z which generates a net magnetic flux H v 
around the vortex. The regular helical vortex strongly 
support the existence of a regular knot. The details will 
be published elsewhere. 

ACKNOWLEDGEMENT 

The work is supported in part by the ABRL Pro- 
gram of Korea Science and Engineering Foundation (R14- 
2003-012-01002-0) and by the BK21 Project of the Min- 
istry of Education. 



[1] P. A. M. Dirac, Proc. Roy. Soc. A113, 60 (1931); Phys. 

Rev. 74, 817 (1948). 
[2] A. Abrikosov, Sov. Phys. JETP 5, 1174 (1957). 
[3] T. H. R. Skyrme, Proc. Roy. Soc. (London) 260, 127 



(1961); 262, 237 (1961); Nucl. Phys. 31, 556 (1962). See 
also, for example, I. Zahed and G. Brown, Phys. Rep. 
142, 1 (1986), and the references therein. 
[4] G. 'tHooft, Nucl. Phys. 79, 276 (1974); A. M. Polyakov, 



9 



JETP Lett. 20, 194 (1974). 
[5] L. Faddeev and A. Niemi, Nature 387, 58 (1997); J. 
Gladikowski and M. Hellmund, Phys. Rev. D56, 5194 
(1997); R. Battye and P. Sutcliffe, Phys. Rev. Lett. 81, 
4798 (1998). 

[6] Y. M. Cho, Phys. Rev. Lett. 87, 252001 (2001). 

[7] Y. M. Cho and D. Maison, Phys. Lett. B391, 360 (1997); 

W. S. Bae and Y. M. Cho, JKPS 46, 791 (2005). 
[8] Y. Yang, Solitons in Field Theory and Nonlinear Analysis 

(Springer) 2001. 
[9] C. Myatt at al, Phys. Rev. Lett. 78, 586 (1997); D. 

Stamper-Kurn, at al, Phys. Rev. Lett. 80, 2027 (1998); 

J. Stenger at al, Nature 396, 345 (1998). 
[10] J. Ruostekoski and J. Anglin, Phys. Rev. Lett. 86, 3934 

(2001); U. Al Khawaja and H. Stoof, Nature 411, 818 

(2001); Phys. Rev. A64, 043612 (2001); H. Stoof at 

al, Phys. Rev. Lett. 87, 120407 (2001); R. Battye, N. 

Cooper, and P. Sutcliffe, Phys. Rev. Lett. 88, 080401 



(2002); C. Savage and J. Ruostekoski, Phys. Rev. Lett. 

91, 010403 (2003). 
[11] J. Nagamatsu et al., Nature 410, 63 (2001); S. L. Bud'ko 

et al., Phys. Rev. Lett. 86, 1877 (2001); C. U. Jung et 

al., Appl. Phys. Lett. 78, 4157(2001). 
[12] H. D. Yang et al., Phys. Rev. Lett. 87, 167003 (2001); J. 

J. Tu et al., Phys. Rev. Lett. 87, 277001 (2001). 
[13] B. Piette, B. Schroers, and W. Zakrzewski, Nucl. Phys. 

439, 205 (1995). 
[14] Y. M. Cho, Phys. Lett. B603, 88 (2004); 

|hep-th/0404T8l1 

[15] Y. M. Cho, |cond-mat/0~ 12325 

[16] Y. M. Cho, |cond-mat/0409636l 

[17] Y. M. Cho, |cond-mat/0112498l |cond-mat /0308182| 
[18] Y. M. Cho, Phys. Lett. B616, 101 (2005). 
[19] Y. M. Cho, |hep-th/0110076| 



